Comment on "Maximal planar networks with large clustering coefficient and 

power-law degree distribution " 
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This Comment corrects the error which appeared in the calculation of the degree distribution of 
random apollonian networks. As a result, the expression of P(k), which gives the probability that 
a randomly selected node has exactly k edges, has the form P(k) oc l/[k(k + l)(fc + 2)]. 

PACS numbers: 89.75.Hc 



In a recent article p|, Zhou et.al. discussed a class 
of networks, called random apollonian networks (RAN), 
which display scale-free degree distribution, very large 
clustering coefficient and very small average path length. 
The simultaneous small-world and scale-free properties 
can mimic well many real-life network systems. 

However the analysis errs in calculating the degree dis- 
tribution function. The rate equation they used to derive 
the expression of P(k) has the form: 



n(N + 1, k + 1) = n(N, k) + n(N, k + 1)(1 
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where n(N, k) denotes the number of nodes with degree 
k when N nodes are present in the network, and the 
total number of triangles. This rate equation describes 
the time evolution of the number of nodes with degree 
k + 1. The first term accounts for the process in which 
a site with k links is connected to the new site, leading 
to a gain in the number of sites with k + 1 links. This 
happens with probability jf—. The second term on the 
right-hand side of Eq. ||TJ accounts for a corresponding 
role (loss). Using the asymptotic limit n(N, k) = NP(k) 
and P(k + 1) — P(k) = dP/dk, they obtained the form 
of the degree distribution function: P(k) oc k 1 with 
7 = (3N + S)/N. However, as we will state below, that 
does not give the right probability of nodes with 
degree k gaining a new link. Thus their result for P(k) 
is incorrect. 

The growth of a RAN is performed by randomly se- 
lecting a triangle in the network, and then a new node 
connects to the three vertices of the triangle (see Ref. [l| 
for details). This means that the correlation among the 
nodes is very strong, i.e., after a new link added to a node, 
the remaining two links are constrained to link to it's two 
neighbors. Thus, on the one hand j^- would include the 
probability of the evolution of other nodes whose degree 
is not equal to k; on the other hand, the evolution of the 
neighbors (except for those neighbors with degree fc) of 
a node (this node has degree k) would effect it's evolu- 
tion in the reverse, which could lead to the information 
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FIG. 1: Schematic illustration of the growth of RAN (left 
column) and AN (right column) for t = 0, 1 and 2 (from top 
to bottom). Each triangle in the RAN is represented by an 
edge in the AN. That a node added to the RAN corresponds 
to an edge is replaced by a triangle in the AN. Left: the red 
node and lines denote the newly added node and links. Right: 
the dash blue line depict the random selected edge (to update) 
and the red lines denote the two newly added edges after each 
generation. Note that the information included in the edges 
are also renewed at the same time (the coordinates of nodes 
in the RAN denoted by the edges in the AN). 



missing for the evolution of n(N, k). The same analysis 
is also suitable for the second term of the Eq. (JIJ. In the 
following, we give out the right form of the expression of 
P(k). 

Instead of n(N,k), let us investigate the time evolu- 
tion of p(k,ti,t), the probability that at time t a node 
i introduced at time U has a degree k. This "node- 
based- method "can overcome the strong node-node cor- 
relation arising in the "group-based-method " . Indeed, 
the growth of the RAN can be understood by the evolu- 
tion of an adjacent network (for convenience, we call it 
AN): each triangle in the RAN is represented by an edge; 
whenever two triangles share a edge, the corresponding 
two edges are connected in the AN; That a new node is 
added to the RAN (leading to two triangles increasing ) 
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N}^ is equal to the degree of the ith node: = fc,-. 

Thus, for ith node, the probability of gaining a new 
link (or the emergence times increasing by 1) is ki/2t, 
where 2t denotes the total edges in the AN at time t (note 
that each updating step would give rise to two extra edges 
increasing) . Consequently the master equation governing 
p(k,ti,t) for the RAN has the form 



r cquE 

a a, 



p(k,ti,t+l) = -^p(k-l,t h t) + (l--)p(k,ti,t) (2) 



with the initial condition p(k, t 
the boundary one p(k,t,t) = Sk,i 
tion can be obtained as Q: 



= 0, 1, t = 1) = Sk,i and 
The degree distribu- 



P(k)= lim (Y.p(k,U,t))/t. 



(3) 



FIG. 2: The distribution of P(k) as a function of A, where 
A=k(k + l)(fc + 2), for different network sizes. The dash line 
is a guide to the eyes with slope -1.0. 



corresponds to a randomly selected edge is replaced by 
a new triangle in the AN (this leads to two edges in- 
creasing) (see figure 1). Note that each edge in the AN 
represents three nodes in the original RAN. When a new 
node enters the system, the degree of node i increases by 
1 with a probability dependent on the emergence times 
of this node on all the edges of the AN, which is equal 
to the number of triangles containing the ith node in the 
RAN. Let denote this number. For large time limit, 



Using Eq. © and the expression p(k)—p(k — l) = dp/dk, 
one can get that P(k) is the solution of the recursive 
equation: P{k) = P(k - l)(fc - \)/{k + 2) for k > m+1, 
and P(m) = 2/(m + 2) for k = m, where m is the degree 
of a node at the time it enters the sytem (in the present 
case m = 3). Solving for P(k) gives that 



P(k) 



2m(m + 1) 



(4) 



fc(fc+l)(fc + 2) 
In Fig. [21 the simulation results of P(k) for differ- 
ent network sizes N = 10 4 , 10 5 and 10 6 are plotted as a 
function of k(k + l)(fc + 2), the best fit line gives out that 
the decay exponent is —1.002 ± 0.003, which is in well 
agreement with the analytic result Eq. 
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